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ABSTRACT 

Context. Polarization is a powerful remote-sensing method to investigate solar system bodies. It is an especially sensitive diagnostic tool to 
reveal physical properties of the bodies whose observational characteristics are governed by small scatterers (dust, regolith surfaces). For these 
objects, at small phase angles, a negative polarization is observed, i.e., the electric vector E oscillates predominantly in the scattering plane, 
contrary to what is typical for rather smooth homogeneous surfaces. The behavior of negative polarization with phase angle depends on the 
size, composition and packing of the scatterers. These characteristics can be unveiled by modelling the light scattering by the dust or regolith 
in terms of the coherent backscattering mechanism. 

Aims. We investigate the surface properties of TNOs and Centaurs by means of polarimetric observations with FORS 1 of the ESO VLT. 
Methods. We have obtained new broadband polarimetric measurements over a range of phase angles for a TNO, 50000 Quaoar (in the R Bessel 
filter), and a Centaur, 2060 Chiron (in the BVR Bessel filters). Simultaneously to the polarimetry, we have obtained R broadband photometry 
for both objects. We have modelled these new observations of Quaoar and Chiron, and revised the modelling of previous observations of the 
TNO 28978 Ixion using an improved value of its geometric albedo. 

Results. TNOs Ixion and Quaoar, and Centaur Chiron show a negative polarization surge. The Centaur Chiron has the deepest polarization 
minimum (-1.5-1.4 %). The two TNOs show differing polarization curves: for Ixion, the negative polarization increases rapidly with phase; for 
Quaoar, the polarization is relatively small (^ -0.6 %), and nearly constant at the observed phase angles. For all three objects, modelling results 
suggest that the surface contains an areal mixture of at least two components with different single-scatterer albedos and photon mean-free paths. 
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1. Introduction 

Transneptunian objects (TNOs) in the Kuiper Belt are consid- 
ered to represent one of the oldest and possibly most original 
population of solar system bodies that can be observed from 
Earth. Centaurs are escapees from the Kuiper Belt through 
gravitational interaction with Neptune and the other giant plan- 
ets. They may eventually become members of the Jupiter fam- 
ily of comets, or may be ejected from the planet region due to 
close encounters with the giant planets. 

The intense study of physical properties of TNOs and 
Centaurs was triggered by the advent of large telescopes on 



* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 69.C-0133 and 073.C-0561 (PI: H. 
Boehnhardt) 



the ground: besides a large set of photometric colours, also 
visible and near-IR spectra of a number of objects are avail- 
able now. Polarimetric observations are more scarce: except 
for Pluto/Charon system (that was observed unresolved, e.g., 
by Kelsey & Fix 1973 1, it was only recently that broadband po- 
larized radiation of a TNO, the Plutino 28978 Ixion, has been 
observed and modelled (Boehnhardt et al. 2004). 

Polarimetry is a powerful tool to investigate the physical 
properties of atmosphereless bodies. At small (< 30°) phase 
angles (the phase angle is the angle between the Sun and the 
observer as seen from the object), these objects exhibit a phe- 
nomenon of negative polarization: the observed flux perpen- 
dicular to the plane Sun-Object-Observer (the scattering plane) 
minus the observed flux perpendicular to that plane, divided by 
the sum of the two fluxes, turns to be a negative quantity. This 
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phenomenon, first discovered through lunar observations by 
Lyot (I1929> . escapes from common sense interpretation, since 
elementary physics tells that reflected electric vector E oscil- 
lates predominantly in the plane perpendicular to the scattering 
plane rather than in the scattering plane. Solar-system objects 
show two types of angular dependence of negative polariza- 
tion: either a smooth phase-angle change that has the minimum 
at ~ 10° (S-, C- asteroids, Moon) or a sharp surge with the 
minimum at ~ 1 — 2° (Saturn rings, Europa, E-asteroids) (see, 
e.g., Rosenbush et al. 2002 1. Both types of negative polariza- 
tion, which also were observed in powdered laboratory sam- 
ples, are currently interpreted in terms of enhanced backscat- 
tering of multiply scattered rays (Shkuratov 1989; Muinonen 
19901. 

Observations of negative polarization and simultaneous 
photometry of main-belt asteroids and other solar system bod- 
ies (see, e.g., Belskaya et al. 2005 , Rosenbush et al. 2005 1 can 
be modelled to infer the properties of the surface texture of 
these objects. Faintness of the targets was the main obstacle 
hampering the same kind of study in TNOs and Centaurs 1 . 
Thanks to the advent of the large telescopes and instruments 
equipped with polarimetric capabilities, observations of TNOs 
and Centaurs are nowadays possible with signal to noise ratio 
comparable to that commonly reached for main-belt asteroid 
observations with small and middle-size telescopes. 

After our first polarimetric study of 28978 Ixion 
(Boehnhardt et al. 12004). in this paper we present new 
polarimetric and photometric measurements obtained with 
FORS1 at the ESO Very Large Telescope (VLT) for a TNO, 
50000 Quaoar, and a Centaur, 2060 Chiron. We also present a 
revised modelling of the observed polarization and photometry 
of 28978 Ixion based on a determination of the geometric 
albedo that has been recently obtained, and that was not 
available at the time of our first modelling effort. 



2. Target summary 

Criteria used for target selection are that the targets are bright 
enough to allow us to measure the polarization with an error 
bar smaller than 0.05 % in less than two hours telescope time. 
With FORS1 at the ESO VLT, this sets the R magnitude limit 
to about 20. Another constraint is the possibility to observe the 
largest possible phase angle range. Complementary informa- 
tion on the geometric albedo and surface composition is essen- 
tial for the modelling part. Moreover, it is desirable to study 
members of the various dynamical and taxonomic groups iden- 
tified among the TNO population (Plutinos, Centaurs, classical 
and scattered-disk Objects). We finally selected three objects 
with known physical parameters (geometric albedos, colours, 
spectral slopes, and surface composition): Ixion, Quaoar, and 
Chiron. 



1 Another difference in the observing and modelling techniques is 
that, due to the larger distance, the observed phase angle range is much 
smaller for TNOs and Centaurs than for main-belt asteroids 



2.1. 28978 Ixion 

28978 Ixion, discovered in 2001, belongs to the dynamical 
class of Plutinos, and it is one of the largest known TNOs (400- 
550 km according to Stansberry et al. 2005|. The visible spec- 
trum by Marchi et al. ( 2003 1 is featureless with a gradient S ' 
of 19.8%/100nm. Optical and near-IR (Licandro et al. 120021 
spectra have been interpreted by Boehnhardt et al. ( 2004) us- 
ing an areal mixture of Titan tholin, amorphous carbon, wa- 
ter ice, and ice tholin. The same authors also present a surface 
model of Ixion based upon their R filter polarimetry and simul- 
taneous R band photometry of the objects spanning the phase 
angle range 0.25° - 1.34°. In that work, an /?-band geometric 
albedo of 0.1 was assumed for the modelling. Here we repeat 
the analysis for the higher R -band geometric albedo now avail- 
able from Spitzer observations (0.23; Stansberry, priv. comm.). 
Further details about the properties of this object are given by 
Boehnhardt et al. d2004b . 

2.2. 50000 Quaoar 

50000 Quaoar is a classical disk object in the Kuiper Belt. 
Orbital elements and red visible colours (Fornasier et al. 2004 1 
suggest that the object could be a member of the "dynamically 
hot" population that is supposed to have migrated to the classi- 
cal disk only after formation closer to the Sun (Gomes 2003 1. 
Apart from Pluto/Charon, 50000 Quaoar is the only TNO so 
far for that disk-resolved photometry could be performed: HST 
measurements allowed to determine the overall size and geo- 
metric albedo of the object to be 1260 ± 190 km, and about 
0.1, respectively (Brown & Trujillo 2004). The photometric 
lightcurve of the object seems to be double-peaked with a pe- 
riod of about 17.6 h and an amplitude of 0.1 3 mag suggesting 
an aspherical shape of the body and/or geometric albedo varia- 
tions of the surface (Ortiz et al. 2003 1. 

Quaoar visible spectra were obtained by Marchi et al. 
( 2003 ) and by Fornasier et al. ( 2004 1. The reflectivity gradients 
S ' obtained in the two papers are not fully consistent, and their 
mean value is 27.6 + 0.3 %/100nm. Visible spectrum appears 
to be featureless. 

Quaoar has been observed also in the near-infrared by 
Jewitt and Luu (2004) at the Subaru 8 m telescope. The com- 
plete spectrum shows a positive-slope continuum from 0.4 up 
to 1 .3 /urn, that is considered typical for the presence of organic 
materials on its surface. The spectrum shows strong absorp- 
tion bands at 1.5 and 2.0yt/m due to H2O ice with the band at 
1 .65 yt/m typical for the crystalline structure in the ice. A small 
presence of ammonia hydrate has also been supposed on the ba- 
sis of the presence of faint features at 2.2 //m (detected also by 
Pinilla-Alonso et al. 2004 1. This was the first time that the qual- 
ity of a TNO spectrum was good enough to distinguish between 
crystalline and amorphous ice. The detection of crystalline ice 
indicates that the temperature has reached at least 110K (crit- 
ical temperature necessary for crystallization). This object is 
large enough to be cry ovolcanic ally active, and crystalline ice 
and ammonia hydrate might be products of this type of activity. 
Jewitt and Luu ( 2004 1 suggested that Quaoar has been recently 



S. Bagnulo et al.: Polarimetry of Kuiper belt TNO objects and Centaurs 



3 



resurfaced either by impacts or by cryovolcanic outgassing or 
by a combination of these two processes. 

2.3. 2060 Chiron 

2060Chiron is the first discovered (1977) and best observed 
Centaur. The intermediate character of the object between 
TNOs and comets is apparent from photometric observations 
that show recurrent episodes of coma activity (gas and dust) 
and of stellar appearance (see for example Meech & Belton 
[T990I Luu & Jewitt 119901 Bus et al. [T99"Tl Duffard et al. l2002l 

Chiron was observed spectroscopically in the visible region 
by many authors (see Barucci et al. 2003 1 showing a flat spec- 
trum with no absorption features. The reflectivity gradient S', 
which ranges between -0.2 up to 2.3 + 0.1 %/100nm, seems 
more similar to that of C-type asteroids rather than to the mean 
reflectance slope of cometary nuclei. The small variation on 
the optical reflectivity gradient could be due to dust production 
variation connected to episodes of recurrent cometary activity. 

Several spectra obtained in the near-infrared did not show 
any features. Only Foster et al. (1999) and Luu et al. (2000) 
detected a 2 pm absorption band suggesting the presence of 
H2O ice on Chiron's surface. Later on, Romon-Martin et al. 



( 2003 1 observed Chiron again during high activity in the visi- 
ble and NIR showing a fiat behaviour without any spectral fea- 
tures finding that is compatible with the hypothesis made by 
Luu et al. (2000 1 that the detection of water ice in Chiron spec- 
tra would be correlated with its cometary activity level. Such 
activity could cause a rain of cometary debris on its surface 
changing the surface mantle. The ice present on the surface, is 
probably mixed with dark impurities which mask the spectral 
bands. 

Nucleus properties using multi-wavelength information has 
revealed an about 70 km nucleus of relatively bright geometric 
albedo (0.17) and moderate axis ratio (1.16) or surface albedo 
variations (Groussin et al. 20041. 

From the large number of publications on Chiron (more 
than 150 to-date) several interesting properties of 2060 Chiron 
have been worked out. However, a synoptic picture of the nu- 
cleus and its surface properties has not yet evolved. 

3. New observations with FORS1 

Observations of Chiron and Quaoar have been obtained at the 
ESO VLT with the FORS1 instrument in service mode dur- 
ing the observing period from April to September 2004. Until 
June 1, 2004, FORS1 was attached at the VLT Unit Telescope 
1 (Antu). After that date, FORS1 was moved to the VLT Unit 
Telescope 2 (Kueyen). 

FORS 1 is a multi-mode instrument for imaging and (multi- 
object) spectroscopy equipped with polarimetric optics. For the 
present study, FORS 1 has been used to measure the broadband 
polarization of Chiron at six different epochs in the Bessel BVR 
filters, and to measure the broadband polarization of Quaoar at 
five different epochs in the Bessel R filter (Sect. 13 . II . In fact, 
one series of Chiron measurements was started on night 2004- 
08-05/06 and aborted after the observations in the R filter be- 
cause the seeing conditions were too good (^ 0.6") with con- 



sequent risk of CCD saturation. The series was repeated the 
following night. From target acquisition images, a by-product 
of polarimetric observations, we could also obtain photometry 
in the Bessel R filter (Sect. l3~2l . Differential tracking was used 
for all our observations, so that long exposure time images of 
Chiron, obtained in polarimetric mode, could be combined al- 
together to search for coma activity (Sect. l3~3l . 

Taking advantage of the flexibility offered by the 'VLT ser- 
vice observing mode', we distributed the observations along 
a few months as to obtain data points approximately equally 
spread over the phase angle ranges of the targets. We set precise 
time intervals for the execution of the observations. In presence 
of the Moon, the sky -background is highly polarized, hence we 
generally tried to avoid observations with the target close to 
the Moon, and with a too large fraction of lunar illumination. 
However, a posteriori we found that the presence of the Moon 
did not jeopardize our observations. The log of the observations 
can be inferred from Tables[2to0] 

3.1. Polarimetry 

To perform linear polarization measurements, a A/2 retarder 
waveplate and a Wollaston prism are inserted in the FORS1 
optical path (see Appenzeller 1 1 967t . The A/2 retarder wave- 
plate can be rotated in 22.5° steps. Stokes Q and U parame- 
ters (defined as in Shurcliff 1962 1 are measured by combin- 
ing the photon counts (background subtracted) of ordinary and 
extra-ordinary beams (/° and / e , respectively) observed at var- 
ious retarder waveplate positions a, where a indicates the an- 
gle between the acceptance axis of the ordinary beam of the 
Wollaston prism and the fast axis of the retarder waveplate. 

In the following, we will always work with the ratios Q/I 
and U/I, and will adopt the notation: 



Q 



and P 



u 



U 

7 



In the ideal case, Pq is obtained measuring the quantity 



r = (-1) 



k r-r 
f°+f e 



at any retarder waveplate position a - k 45°, and Py is ob- 
tained measuring the ratio r at any position a = k45° + 22.5° 
(k = 0, 1, 2, . . . , 7). The validity of this assertion can be veri- 
fied e.g. with the help of Eq. (1.33) of Landi DeglTnnocenti & 



Landolfi (2004i. In practice, there are several deviations from 
the ideal case. For instance, the actual retardance value of the 
retarder waveplate may deviate from the nominal n value; the 
transmission of the ordinary and extraordinary beam are not 
identical, even after flat fielding correction. The effect of these 
(and other) sources of instrumental polarization can be largely 
reduced at the first order by measuring 



p'j 

v 



(£=£) _(f^Jl) 

V/°+/ e /a=2(i-l)x45° \f°+f°/a-- 



(2j-l)x45° 



\f°+f°Ja=2(i 



2(i-l)x45°+22.5° 



(1) 



(2;-l)x45°+22.5° 
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where i and j are integers numbers 2 . In the simplest case, linear 
polarization can be measured from the observations obtained at 
four angles of the retarder waveplate: 



Pn = 



(/"+/* )o=0 o (f°S )«=45° j 
\/ u +/ c Ar=22.5° \f°+f c )a=67.5 



It is convenient (and recommended in the FORS1/2 user 
manual) to obtain Stokes Q and U adding up observations ob- 
tained with the retarder waveplate at various positions: 



p _ 1 v m p" 
- m 2w=i r Q 

p _ 1 v m p" 



(2) 



where m represents the number of pairs of observations for each 
Stokes parameter, and Pq and P l v are obtained from Eq. 
setting i = j — I. We performed simple numerical simulations 
to study the impact on the precision of the polarimetric mea- 
surements of a deviation of the waveplate retardance from its 
nominal value (180°). We found that using Eq. (0 with m = 2, 
a deviation from the nominal value of the waveplate retardance 
as large as 5°, for the polarization value observed in our targets, 
would introduce a spurious contribution <K 0.01 %. Figure 4.1 
of FORS1/2 user manual shows that actual deviation of the re- 
tarder waveplate from the nominal value is well within 5°. We 
conclude that in our data, the effect of instrumental polarization 
due to the chromathism of the retarder waveplate is negligible. 

It should be noted that FORS1 is affected by a problem 
with spurious linear instrumental polarization that cannot be 
eliminated even by using the reduction technique explained 
above. This spurious polarization, due to the presence of rather 
curved lenses in the collimator, combined with the non com- 
plete removal of reflections by the coatings, is axially symmet- 
ric, and smoothly increases from less than 0.03 % on the optical 
axis to 0.7 % at an axis distance of 3 arcmin (in the V band). 
This problem has been discovered and investigated by Patat & 
Romaniello (2005 1. In our case, since our targets are always in 
the center of the field of view, the problem of the instrumental 
polarization can be safely ignored. 

The error-bar due to photon-noise on the Pq or Py mea- 
sured from a pair of observations is 



((^) 2 4 + ((7fe) 2 4L, 
((^) 2 4 + ((7fe) 2 4L (2/ 



where <f>o = in case X represents Pq and <p$ = 22.5° in case 
X represents Py. When Pq and Pu are obtained adding up m 
pairs as in Eq. 0, the error is given by 



1 m 

— V a 2 ,, 
m j-J x 



2 We note that in a similar formula that we reported in Boehnhardt et 
al. |2004 1, the factor 2 in the indices that denote the angles is missing 
because of a typo. 



For the sake of simplicity, assuming f°-f e -N and also as- 
suming N to be the same for all positions of the retarder wave- 
plate, we obtain 



1 1 



(4) 



where V77 represents the signal to noise ratio (SNR) measured 
in the individual beam for each of the m exposures (in other 
words, yfmN is the cumulative SNR in each beam). It appears 
for instance that to get an error bar on the Stokes parameter Pq 
or Pu of about 0. 1 %, one should take a pair of exposures with 
SNR of about 500 each (integrated over the individual point- 
spread function area of each beam of each exposure). 

When multiple pairs of exposures are taken, it is useful to 
study the distribution of the P'l (P' v ) values obtained substitut- 
ing i,j = 1,2, . . .m. In particular, a cr-clipping algorithm can 
be applied to the P l L and P' v distributions in order to "clean" 
the data, by rejecting those values that deviate more than a cer- 



tain distance from the median (see Boehnhardt et al. 2004 and 
Bagnulo et al. 120051 . 

Stokes Q and U are usually measured with the instrument 
position angle = 0°, i.e. to have the acceptance axis of the ordi- 
nary beam of the Wollaston prism aligned to the North Celestial 
Meridian (and the acceptance axis of the extra-ordinary beam 
perpendicular to it). We then transform the Stokes parameters 
according to 



P' Q = cos(2(<£ + n/l)) P Q + sin(2(<D + tt/2)) P v 
P'jj = - sin(2(<D + tt/2)) Pq + cos(2(<D + tt/2)) P v 



(5) 



where O is the angle between the direction Object-North Pole 
and the direction Object-Sun. This angle can be calculated ap- 
plying the four parts formula to the spherical triangle defined 
by the object (with coordinates aT,<5x), the Sun (with coordi- 
nates (ar ,£ ) and the North celestial pole: 



sin <St cos(a - aj) = cos(<St) tan(5 Q ) - sin(a - aj) 



1 



tan(O) 



This way P'g represents the flux perpendicular to the plane Sun- 
Object-Earth (the scattering plane) minus the flux parallel to 
that plane, divided by the sum of the two fluxes. The angle of 
maximum polarization is obtained as 



©' = - arctan 
2 



(P'r 



p: 



(6) 



Q> 



y-') where 

-l)x45°+</. 



or 



0' 



if P' Q > and P' v > 
if P' e >0andP' £/ <0 
ifPg<0 



Jtt/4 if P'g = OandP'^ > 
| 3tt/4 ifPg ^OandP'y < 



(see Landi DeglTnnocenti & Landolfi 2004 ) Incidentally, it 
should be noted that the 0^ term is occasionally incorrectly 
neglected. 
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Observations of Chiron were performed with the retarder 
waveplate at all positions between and 157.5° (at 22.5° steps), 
i.e., setting m — 2 in Eq. 0, using the broadband Bessel B, V, 
and R filters. For Quaoar we used all positions of the retarder 
waveplate from 0° to 337.5°, i.e., setting m — 4 in Eq. (|2j, us- 
ing the Bessel R filter. For each frame, exposure times t were 
as follows. For Chiron we set t = 460, 140, and 110 s in the 
B, V, and R filter, respectively. For Quaoar (R filter only), we 
set t = 250 s. Note that, for each Stokes parameter we obtained 
2 pairs of exposures for Chiron, and 4 pairs of exposures for 
Quaoar. Therefore the total shutter time for Chiron observa- 
tions was of about 60, 20, and 15 min in the B, V, and R filter, 
respectively, and of about 67 min for Quaoar (R filter only). 
For each observation of Chiron, the typical SNR accumulated 
in each beam was about 1700, 1600, and 1600, in the B, V, and 
R filters respectively. For each observation of Quaoar, in the R 
filter, we obtained a SNR accumulated in each beam of about 
1400. 

The photon counts f° and f e were measured via simple 
aperture photometry performed on the images, obtained after 
bias subtraction and flatfield correction (master fiat field was 
obtained from sky images obtained during twilight with no 
polarimetric optics in). More sophisticated methods based on 
point spread function (PSF) fitting are difficult to apply because 
of the star trailing due to differential tracking on the moving 
targets. Sky background was measured in an annulus with 5 " 
inner radius, centered around the source, of 2" to 6" width. 
The errors on f° and / e were estimated as explained in the 
Sect. 3.3.5.8 of Davis ( 1987). f°, f e , and their error estimates 
were measured for aperture values ranging from 0.6"to 4". 

Pq and Pjj values were found to be slightly dependent on 
the aperture adopted to measure f° and f e . This effect was 
more critical when the target was in a crowded field and/or in 
presence of strong background polarization by the Moon. The 
final value was selected as the one for which the error on Pq 
and Pjj was minimum, i.e., usually for aperture values between 
1.2" and 2.5". 

The results of the polarimetric measurements are given in 
Tables^andlJ] that are organized as follows. Columns 1 and 2 
give the epoch of the observations (date and UT time). 

Columns 3 and 4 give Moon angular distance and fraction 
of lunar illumination (FLI), respectively. 

Column 5 gives an estimate of the night time sky condi- 
tions: THN — thin cirrus, CLR= clear, PHO = photometric. The 
classification given in these Tables is based on our inspection 
of reduction products and of the atmospheric monitors of the 
observatory at 

http : //www . eso . org/gen-f ac/pubs/astclim/ 
forecast/meteo/CIRA/images/repository/lossam/ 
and is somewhat arbitrary. Those nights when no photomet- 
ric standard stars were observed, or when the zeropoints were 
not considered stable for the entire night, are indicated with a 
question mark. Sky transparency does not affect the precision 
of the polarization measurements, but it does affect the preci- 
sion of the photometry (see Tables[3]and|4j. Note that the pre- 
cision of the measurements (both polarimetric and photomet- 
ric) depends also on how the field of view close to the target 
is crowded with background objects. This situation of course 



changes from epoch to epoch. Observations labelled with (c) 
are hampered to some extent by a crowded background. In the 
case of Chiron observations on 2004-09-01 the situation was 
probably especially critical. 

Column 6 gives the Sun-Target-Observer angle, i.e., the tar- 
get's apparent phase angle as seen at observer's location, ex- 
pressed in degrees. This phase angle was obtained from the 
object ephemeris calculated at the JPL's solar system dynamics 
WWW site at http://ssd.jpl-nasa.gov 

Column 7 specifies the broadband filter used for the obser- 
vations. 

Columns 8, 9, and 10 give Pq , Pjf , and the angle of max- 
imum polarization 0' , respectively, after the transformation of 
Eq. (jSJi. It should be noted that for both Chiron and Quaoar, the 
measured PC, is generally consistent with zero (equivalent to 
the fact that 0' is generally consistent with 90°). This means 
that the principal axes of the polarization ellipse are aligned 
with the coordinate axes perpendicular and parallel to the scat- 
tering plane. Furthermore, the observed P'q, i.e., the flux per- 
pendicular to the scattering plane minus the flux parallel to that 
plane, divided by the sum of the two fluxes, is always a nega- 
tive quantity. This means that the direction of the polarization is 
included in the scattering plane (this case is normally referred 
to as "negative polarization"), in contrast to what is expected 
for a dielectric medium, i.e., that the direction of polarization 
be perpendicular to the scattering plane (this case is normally 
referred to as "positive polarization"). 

Polarimetric measurements of Quaoar and Chiron are plot- 
ted against object phase angle in Figs.[3]and|4] As far as Chiron 
observations are concerned, it appears that the observed polar- 
ization does not depend on the filter used, at least when consid- 
ering the error bars of the observations. The absolute value of 
the polarization increases as the phase angle decreases, perhaps 
reaching a minimum around phase 1.5-2.0 deg. Observations 
of this object at smaller phase angles were unfortunately not 
taken although originally scheduled in our service observing 
campaign. 

3.2. Photometry 

The object magnitudes were measured in the acquisition im- 
ages obtained in the R Bessel filter with no polarimetric optics 
in the light path. Exposure times were 5 s for Chiron, and 30 s 
for Quaoar. 

Since our program was aimed at polarization measure- 
ments, we did not obtain a number of observations of photo- 
metric standard stars sufficient to estimate precise zeropoints 
and extinction coefficients for the various nights. In most of the 
cases, the zeropoints were obtained from only one frame ob- 
tained during the night when our observations were executed 
(and not necessarily close to them in time). In some cases, no 
photometric standard stars were observed at all (see Tables ^ 
and |3- In these cases we adopted the values measured in the 
nearest night with similar sky conditions. For the extinction 
coefficient Kr and the colour term Icvr we used the values esti- 
mated for the P73 ESO observing period published at 
http : //www . eso . org/observing/dfo/quality/FORSl/ 
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Table 1. The observed polarization of 2060 Chiron. The meaning of the various columns is given in the text 



Date Time(UT) Moon Sky Phase angle Filter P' Q F v & 



(yyyy mm dd) 


(hh:mm) 


dist. 


FLI 


transp. 


(DEG) 




(%) 






(%) 




(DEG) 




2004 05 27 


06: 


:37 


139° 


0.3 


CLR 


3 


.469 


R 


-1.14 d 


= 


.04 


0. 


.02 = 


zO 


.04 


89 


.4 = 


z 1 


.0 


2004 05 27 


07 


:03 








3 


.468 


V 


-1.20d 


z 0. 


.04 


0. 


.13 = 


zO 


.04 


86 


.9 = 


zO 


.9 


2004 05 27 


07 


:52 








3 


.467 


B 


-1.24 d 


zO. 


04 


-0. 


.04 = 


zO 


.05 


90 


.9 = 


z 1 


.1 


2004 06 29 


03 


:19 


64° 


0.4 


CLR? 


1 


.410 


R 


-1.40 = 


= 0. 


.05 


0. 


.05 = 


zO 


.04 


89 


.1 = 


zO 


.8 


2004 06 29 


03 


:44 








1 


.409 


V 


-1.35 d 


zO. 


.05 


-0. 


.08 = 


zO 


.04 


91 


.6 = 


zO 


.9 


2004 06 29 


04 


:33 








1 


.407 


B 


-1.32 d 


zO 


.04 


0. 


08 = 


zO 


.05 


88 


.2 = 


z 1 


.2 


2004 08 06 


05 


:53 


94= 


0.7 


PHO 


1 


.766 


R 


-1.40 = 


zO. 


.04 


0. 


.00 = 


zO 


.04 


89 


.9 = 


zO 


.8 


ZUU4 Uo U / 


\)5 


.40 


tlD 


U, / 


TTT\T 
IriiN {C) 


1 
1 


.0/6 


D 

K 


— 1.33 = 


z U. 


.U4 


— U. 


.Ul d 


z U 


,U4 




,1 j 


z U 


Q 
.O 


1(\C\A HQ CM 
ZUU4 Uo U / 


04 


:13 








1 


.829 


V 


-1.46 = 


zO 


A Li 


— u. 


14 d 


z 




92 


.7 -: 


z 


7 


2004 08 07 


05 


:00 








1 


.832 


B 


-1.52 = 


zO 


.04 


-0. 


.09 d 


zO 


.05 


91 


.8d 


zO 


.9 


nnn/l HQ. 1 T 
ZUU4 Uo 1.5 




■ 1 7 


to / 


u.y 


pun /^,\ 
rriU {C) 


L 


O 1 Q 

.ZLy 


r> 
K 


i ii j 
— l.jj = 


z U. 


C\A 
.U4 


A 

u. 


.Uo = 


A 
z U 


n i 




i . 

.5 = 


z 1 


.U 


nnn/1 HQ 1 T 
ZUU4 Uo 1.5 


00 


:42 








2 


.221 


1/ 
V 


-1.34 = 


z 0. 


.U4 


— U. 


.08 = 


z 


C\A 
,U4 


91 


.6 z 


z 


o 


O AA a no 1 1 

2004 08 13 


01 


:31 








2 


.223 


B 


-1.45 = 


zO. 


.03 


-0. 


.19 = 


zO 


.05 


93 


.8 = 


z 1 


.1 


2004 09 01 


00 


:56 


74° 


0.5 


THN? (c) 


3 


.327 


R 


-1.12 = 


zO 


.04 


0. 


.05 = 


zO 


.04 


88 


.7 = 


z 1 


.1 


2004 09 01 


01 


:21 








3 


.328 


V 


-1.09 = 


zO 


.05 


-0 


08 d 


zO 


.04 


92 


.1 d 


z 1 


.2 


2004 09 01 


02 


:10 








3 


.330 


B 


-1.15 = 


zO. 


.05 


-0. 


.24 = 


zO 


.06 


95 


.9 = 


= 1 


.6 


2004 09 27 


01 


:24 


56° 


0.4 


CLR? 


4 


.232 


R 


-0.89 = 


zO. 


.07 


0. 


.03 = 


zO 


.05 


89 


.1 d 


= 1 


.8 


2004 09 27 


01 


:49 








4 


.232 


V 


-1.05 = 


zO 


.05 


-0. 


.20 = 


zO 


.05 


95 


.4 = 


z 1 


.7 


2004 09 27 


02 


:38 








4 


.233 


B 


-0.87 = 


zO 


.09 


-0. 


.10 = 


zO 


.08 


93 


.3 = 


= 3 


.0 



Table 2. The observed polarization of 50000 Quaoar. The meaning of the various columns is given in the text 



Date 


Time (UT) 


Moon 




Sky 


Phase angle 


Filter 


P' 


P'u 


0' 


(yyyy mm dd) 


(hh:mm) 


dist. 


FLI 


transp. 


(DEG) 




(%) 


(%) 


(DEG) 


2004 04 18 


04:43 


47° 


1.0 


CLR 


0.952 


R 


-0.64 ± 0.05 


0.14 ±0.08 


83.9 ±3.9 


2004 05 13 


07:38 


93° 


0.8 


CLR? 


0.496 


R 


-0.50 ± 0.06 


-0.14 ±0.05 


97.7 ± 3.6 


2004 05 26 


03:37 


107° 


0.2 


PHO 


0.252 


R 


-0.49 ± 0.06 


0.10 ±0.05 


84.3 ± 3.5 


2004 07 10 


01:44 


118° 


0.7 


PHO (c) 


0.797 


R 


-0.53 ± 0.04 


0.09 ± 0.05 


85.3 ±3.1 


2004 08 1 1 


02:34 


75° 


0.8 


CLR? 


1.231 


R 


-0.65 ± 0.04 


-0.01 ±0.06 


90.4 ± 2.6 



qc/photcoef f /photcoef £s_for s 1 . html 
i.e., 

K R = 0.045 + 0.019 
k VR = 0.0090 ±0.0017 

For the Quaoar colour indices we adopted V - R = 0.64 ± 0.01 
(Tegleretal. l2003l . and for Chiron V-R = 0.37 ±0.03 (Davies 
et al. 1998 1. Finally, reduced magnitude R of the target was 
obtained by taking into account the distance Sun-Object r and 
the distance Earth-Object A at the epoch of the observations, 
and calculating the reduced magnitude 

R = m R - 5 log(r) - 5 log(A) (7) 

where m R is the observed /?-band magnitude. 

The results of our photometry are given in Tables[3]and|4] 
From a linear least-square fit of the data, we found that 
for Chiron the Hr absolute magnitude at zero phase angle 
is 5.52 + 0.07 mag, and the slope of the opposition surge is 
0.045 + 0.023 mag/deg. Previous numerous photometric obser- 
vations of Chiron (for summary see Groussin et al. 2004 1 have 
shown considerable variations of absolute magnitude with he- 
liocentric distance attributed to its sporadic cometary behav- 



ior. The absolute magnitude of Chiron in 2004 is brighter as 
compared to the observations in 1988-1990 made at the same 
heliocentric distances (Bus et al. I199U . It may indicate that 
we observed Chiron closely after an activity period which may 
have caused partial resurfacing of the object. This is why we 
paid a special attention to search for a coma around Chiron 
(see Sect.l3~3l. 

For Quaoar, the Hr absolute magnitude at zero phase angle, 
obtained from a linear fit is 2.16 + 0.05 mag, and the linear 
slope is 0.16 + 0.06 mag/deg. The determined value of absolute 
magnitude coincides with that used for Quaoar for geometric 
albedo determination (Brown & Truiillo l2004> . 

Note that Hr absolute magnitudes at zero phase angles and 
linear slopes will be re-estimated based on our modelling tech- 
nique in Sect.0] 

3.3. Search for a coma around 2060 Chiron 

Visual inspection and measurements of the full-width-at-half- 
maximum of the Chiron and neighbouring star images did not 
reveal any indications for the presence of a coma around the 
Centaur. For a more thorough analysis we calculated the radial 
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Table 3. The observed R photometry of 2060 Chiron . — , — , — , — , — , — , — , — , — n 



Date 


Time (UT) 


Phase angle 




R 


(WW mm dd) 


(hh:mm) 


(DEG) 


(mag) 


(mag) 


2004 05 27 


06:23 


3.470 


16.64 ± 0.04 


5.68 


2004 06 29 


03:04 


1.411 


16.48 ± 0.05 


5.56 


2004 08 06 


05:40 


1.766 


16.49 ± 0.04 


5.54 


2004 08 07 


03:33 


1.828 


16.65 ± 0.04 


5.69 


2004 08 13 


00:02 


2.219 


16.55 ± 0.05 


5.58 


2004 09 01 


00:39 


3.326 


16.65 ± 0.04 


5.63 


2004 09 27 


01:10 


4.232 


16.84 ± 0.04 


5.74 


Table 4. The observed R photometry of 50000 Quaoar 


Date 


Time (UT) 


Phase angle 


m R 


R 


(yyyy mm dd) 


(hh:mm) 


(DEG) 


(mag) 


(mag) 


2004 04 18 


03:59 


0.952 


18.60 ±0.05 


2.27 


2004 05 13 


06:54 


0.497 


18.60 ±0.04 


2.27 


2004 05 26 


02:52 


0.252 


18.51 ±0.04 


2.19 


2004 07 10 


00:59 


0.796 


18.57 ± 0.04 


2.24 


2004 08 1 1 


01:49 


1.231 


18.73 ± 0.04 


2.38 



27/05/2004 



29/06/2004 



27/09/2004 




Fig. 1. Coma search around 2060 Chiron. Solid lines show the 
Chiron radial profiles, the dotted lines show the radial profiles 
of the comparison stars. The x:— axis represents the logarithmic 
distance in pixels from the centre of the Chiron (or the stars) 
PSF images. The y-axis represents the flux (in arbitrary units). 



Table 5. The best fit coherent-backscattering model param- 
eters for Ixion, Quaoar, and Chiron. We give the single- 
scattering albedos u> and dimensionless mean free paths k( for 
the dark (subscript d) and bright components (b), the weight 
of the dark component wj, the rms values of the polarimetric 
fits, as well as the 7?-band absolute magnitudes Hr and slope 
parameters Icr. 



28798 Ixion R Bessel polarimetry 

-i i i i r 





Ixion 


Quaoar 


Chiron 




0.45 


0.35 


0.15 




250 


300 


120 




0.80 


0.50 


0.60 


4 


20 


10 


500 


Wd 


0.74 


0.46 


0.14 


rms 


0.029 % 


0.069 % 


0.067 % 


Mr 


3.25 


2.15 


5.41 


k R 


0.12 deg" 1 


0.11 deg" 1 


0.041 deg" 1 



flux profile measured in concentric rings around the object, and 
compared it with the radial profile across the star trails, aver- 
aged along the trail direction and for both sides of the trail. 
The resulting object and star profiles are normalized to unity 
brightness at 'one pixel' central distance and to zero at back- 
ground distance 30 pixels from the center. Figure ^shows the 
results for three observing dates (for other dates the presence of 
star trails close to the object image jeopardized the accuracy of 
this analysis). Since the radial profile of Chiron (solid line) is 
basically identical to that of the comparison stars (dotted line) 
or falls at slightly lower flux levels, we conclude that a coma 
around Chiron is not present or it is well beyond our detection 
limit (order of 3 1 mag/arcsec) - if present at all. Hence, we as- 
sume that the polarimetry of the object is not contaminated by 
dust around the object. 



3- 



-0.5 - 




-1.5 - 



PHASE ANGLE (DEG) 

Fig. 2. The observed broadband polarization of 28978 Ixion. 
The solid symbols refer to PL, i.e., the flux perpendicular to the 
plane Sun-Object-Earth minus the flux parallel to that plane, 
divided by the sum of the two fluxes. The empty circles refer to 
PL (measured in the scattering plane as explained in the text). 
Solid line represents the modelling fit to PL. The dotted line 
shows the expected P' u null values. 



8 



S. Bagnulo et al.: Polarimetry of Kuiper belt TNO objects and Centaurs 



50000 Quaoar R Bessel polarimetry 




1 2 

PHASE (DEG) 



Fig. 3. The broadband polarization of 50000 Quaoar observed 
in the Bessel R filter. As in Fig.|2] solid symbols refer to PL, and 
open symbols refer to P' v . Solid line represents the modelling 
fit to P'q. Dotted lines show the expected null P' v values. 

4. Modelling of the observed polarization 

We interpret the polarimetric and photometric phase curves of 
Ixion, Quaoar, and Chiron through extensive numerical simu- 
lations of coherent backscattering by Rayleigh scatterers. Note 
that we avoid making an assumption that the fundamental scat- 
terers responsible for the coherent backscattering contribution 
would be the single particles in the regolith. With the phe- 
nomenological modeling currently including the first multipole 
contribution of an electric dipole, we allow for the possibil- 
ity that the fundamental scatterers can be the volume and sur- 
face inhomogeneities within and on the particles, respectively. 
Shadowing among the regolith particles is known to contribute 
to the opposition effect but not to the negative polarization 
surge. Here, as in Boehnhardt et al. (2004 1, the shadowing ef- 
fect manifests itself in the residual slope parameter resulting 
from the combined coherent backscattering modeling on the 
polarimetry and photometry. 

The coherent-backscattering mechanism is a multiple- 
scattering mechanism for scattering orders higher than the first 
one. For a recent review of the coherent-backscattering mech- 
anism and its relevance in asteroid studies, see Muinonen et 
al. (|2002i and Muinonen (2004i. Note that the computational 
technique for coherent backscattering accounts for all orders of 
scattering contributing in a non-negligible way to the backscat- 
tering peaks and polarization surges. 



2060 Chiron BVR Bessel polarimetry 

i i i i 




1 2 3 4 5 

PHASE ANGLE (DEG) 

Fig. 4. The observed broadband polarization of 2060 Chiron. 
As in Fig. 13 the solid symbols refer to P'q, and the empty sym- 
bols refer to PL. Triangles, squares, and circles refer to the po- 
larization measured in the Bessel B, V, and R filters, respec- 
tively. Solid line represents the modelling fit to P' Q for the R 
filter. The dotted line shows the expected P' v null values. 



In the following we limit ourselves to illustrate the mecha- 
nism at the second order of scattering. Let us consider a semi- 
infinite random medium of discrete scatterers, constrained by 
a plane-parallel boundary with free space, and an electromag- 
netic plane wave incident on the random medium from the free 
space. Let us assume that the incident wave interacts with one 
of the scatterers, giving rise to first-order scattering, and that, 
subsequently, the first-order scattered field interacts with an- 
other scatterer, giving rise to second-order scattering. Let us 
assume that the field scattered by the second scatterer escapes 
the medium and is detected by the observer in the free space. 
Now there is the reciprocal sequence of scatterings in the op- 
posite direction involving the very same two scatterers. In the 
exact backward scattering direction, due to the identical lengths 
of the propagation paths, the two reciprocal wave components 
interfere constructively whereas, in other directions, the inter- 
ference characteristics vary. After configurational averaging, a 
backscattering peak results in the proximity of the backward 
scattering direction. Focusing in on the polarization character- 
istics, the interference is selective and tends to invert the linear 
polarization characteristics in single scattering: for positively 
polarizing single scattering, coherent backscattering tends to 
result in negative polarization. The illustration can be readily 
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3.1 



28978 Ixion R Bessel photometry 



50000 Quaoar R Bessel photometry 



3.2 - 




2.1 - 




PHASE ANGLE (DEG) 



PHASE ANGLE (DEG) 



Fig. 5. The observed reduced magnitude R of Ixion. The 
dashed line represents the best fit obtained with a straight line. 
The solid line represents the model fit. The scattering around 
the fits is fully within the amplitude of intrinsic variability of 
the object. 

generalized to higher orders of scattering. In the modeling that 
follows, all relevant orders of scattering are taken into account. 

For the modelling of Ixion, Quaoar, and Chiron, we car- 
ried out coherent backscattering computations for a total of 360 
spherical media of Rayleigh scatterers, as follows: we assumed 
18 different single-scattering albedos 

cj = 0.05, 0.10,..., 0.90 

and 20 different dimensionless mean free paths 

kl = 2n£/A = 10, 20, 30, ... , 100, 120, 140, 
200,250,300, 400,500 

where k and A are the wave number and wavelength, respec- 
tively. 

As to the /?-band geometric albedos we adopted pa = 
0.23 for Ixion (Stansberry, priv. comm.), 0. 1 1 for Quaoar (cal- 
culated using the diameter from Brown & Trujillo 2004 and our 
determination of the absolute magnitude), and 0.17 for Chiron 
(Barucci et al. I2004> . In our previous interpretation of Ixion 
polarimetric data (Boehnhardt et al. 2004 1, we had assumed 

PR = 0.1. 

We compared the polarimetric observations against the 
spherical media composed of monodisperse Rayleigh scatter- 
ers with given single-scattering albedo and mean free path (two 



Fig. 6. The observed reduced magnitude R of Quaoar. Symbols 
are the same as Fig.|5] 



parameters). The fits were poor: for a fixed geometric albedo, 
the monodisperse Rayleigh-scattering model tends to result in 
polarizations that are substantially pronounced as compared to 
the polarizations observed. 

As in Boehnhardt et al. (2004i, we then studied a two- 
component Rayleigh-scattering model consisting of dark and 
bright scatterers. There are five parameters in such a model: 
two single-scattering albedos and two mean free paths, and the 
weight factor for the dark component (one minus the weight 
factor of the bright component). Fixing the geometric albedo 
fixes the weight factor, reducing the number of free parameters 
to four. After a systematic study of physically realistic combi- 
nations of the two kinds of scatterers, satisfactory fits were ob- 
tained for the polarizations of all three objects. For Ixion, the 
data point at phase angle 0.43° was omitted as an outlier. The 
model parameters and rms values of the fits are summarized in 
Tableland the actual fits are depicted in Figs.[2]-|4] After the 
polarization fits, approximate brightness fits were obtained by 
varying the absolute magnitude Hr and slope parameter of a 
linear phase dependence multiplying the coherent backscatter- 
ing contribution (see Fig.|5]-0. For Chiron, the observation at 
1.8° phase angle was omitted as outlier. 

For all three objects, the dark component shows a mean 
free path substantially longer than the wavelength. For Chiron, 
the mean free path of the bright component is also consider- 
ably longer than the wavelength. The dark components of Ixion 
and Quaoar resemble one another and the bright component of 
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5.4 Mm 




5.8 - 



5 ' 9 1 2 3 4 5 

PHASE ANGLE (DEG) 

Fig. 7. The observed reduced magnitude R of Chiron. Symbols 
are the same as Fig. [5] 

Chiron: this is concluded from the similarity of the parame- 
ters of the corresponding scatterers on the surfaces of the three 
objects. The phase curves of Ixion and Chiron resemble each 
other to an extent where their combined polarimetric and pho- 
tometric data could be explained using a single two-component 
scattering model. 

The polarimetric observations suggest that Ixion and 
Chiron could have similar surface structure. It is notable that 
the geometric albedo ranges of Ixion and Chiron overlap, 
whereas Quaoar stands out as a darker object. More observa- 
tional and theoretical work is required for further conclusions. 

5. Discussion and Conclusions 

We have presented the results of the first polarimetric ob- 
servations for a Centaur and a classical disk object in the 
Kuiper Belt. Together with polarimetric observations of Ixion 
(Boehnhardt et al., 2004 1, a representative of the Plutino pop- 
ulation in the belt, they give us a first idea about the behav- 
ior of linear polarization and intensity of the light reflected 
by surfaces of Kuiper belt objects and Centaurs. Two Kuiper 
Belt objects, Ixion and Quaoar, were observed practically in 
the same range of phase angles (0.25°-1.3°). They have shown 
completely different polarization-phase behaviour. For Ixion 
the negative polarization rapidly increases with phase angle 
(- 1 .02 ±0.025 %/deg), while for Quaoar the polarization degree 
is changed very slowly (-0.17 + 0.10 %/deg). Our polarimet- 
ric observations of Chiron were made for larger phase angles, 



1.4°-4.2°. They are characterized by a pronounced branch of 
negative polarization with a minimum of 1 .4-1 .5 % at phase an- 
gles of 1.5°-2.0°. Such polarimetric characteristics are unique 
among Solar System bodies. For the majority of Solar System 
bodies the polarization at the phase angle 1° is characterized 
by values within 0.1-0.5%. The largest values of the polar- 
ization at 1° measured for non-TNO objects were 0.83 % for 
the dark side of Iapetus (albedo = 0.05) and 0.73 % for Saturn 
Ring A (albedo = 0.75) (Rosenbush et al. 120021 that is no- 
ticeably smaller than the values obtained at our observations 
of Ixion and Chiron. Polarimetric observations of Chiron at 
smaller phase angles are urgently needed for a better mod- 
elling, and to compare Chiron data with the polarimetric curve 
of Ixion and Quaoar. Further observations would be also desir- 
able to identify the inversion angle. 

Our modelling has shown that the possible way to explain 
observed polarization properties of KBOs and Centaur is to 
assume two-component surface media consisting of dark and 
bright scatterers. Such a model succeeds in fitting all observed 
polarimetric and photometric characteristics of the three ob- 
jects and the two-component modelling is realistic. A more 
thorough theoretical study is beyond the scope of the present 
article and should be carried out in the nearest future. 
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help to identify adequate observing periods for our targets (i.e., with 
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